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ABSTRACT 

We explore the impacts of ultraviolet (UV) radiation feedback on galaxies during 
the epoch of reionisation by cosmological simulations in which hydrodynamics and 
the transfer of the H and He ionising photons are consistently coupled. Moreover we 
take into account H2 non-equilibrium chemistry, including photo-dissociation. The 
most striking feature of the simulations is a high spatial resolution for the radiative 
transfer (RT) calculation which enables us to start considering not only external UV 
feedback processes but also internal UV feedback processes in each galaxy. We find that 
the star formation is significantly suppressed due to the internal UV and supernova 
(SN) feedback. In low mass galaxies with M < 1O 9 M , a large amount of gas is 
evacuated by photo-evaporation as previous studies have shown, which results in the 
suppression of star formation. Surprisingly, star formation in massive halos is also 
strongly suppressed despite the fact that these halos hardly lose any gas by photo- 
evaporation. The suppression of star formation in massive halos is mainly caused by 
following two factors; (i) small scale clumpy structures in the galaxies are smoothened 
by the internal feedback, (ii) although the dense gas in the galaxies is mostly neutral, 
the H2 formation and cooling processes are disturbed by mild photo-heating. Photo- 
dissociating radiations actually suppress star formation, but the magnitude of the 
effect is not so large in massive galaxies. Even though our simulation volume is too 
small to be a representative patch of the Universe during reionisation, we find that 
our simulated star formation rate densities and HI fractions at z ~ 6 — 7 are consistent 
with those found in observations. 
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1 INTRODUCTION 

Many observations have indicated that the Universe was 
reionised early on. Acc ording to the Gunn- Peterson effect 
(|Gunn &; Petersonl 1 19651 ) in the spectra of high-; ? quasars, 
the U niverse was already reionised by z ~ 6 (|Fan et al.l 
2006a. b). In addition, the Wilkinson Microwave Anisotropic 
Probe (WMAP) satellite produced further information on 
the reionisation epoch. The measured optical depth for the 
Thomson scattering of CMB photons by free electrons is 
r e = 0.087 =b 0.014. Assuming an instantaneous reionisation 
as a toy model, this value indicates that the Universe was 
reionised at z = 10. 4 ±1.2 ([jarosik et al.ll201lh . Moreover, in 
the last decade, some observed quantities, such as the spec- 
tra of Gamma-ray burst 050904 (|Totani et alj f2006) and the 
evolution of the luminosity function of Lyman-a emitters 
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(LAEs) (|Ouchi et al.ll2009h . have also been used to evaluate 
the neutral hydrogen fraction during the epoch of reionisa- 
tion. However the exact reionisation history is still unclear, 
since we cannot directly detect the ionisation state of inter- 
galactic medium (IGM) even with the latest observational 
facilities. 

Numerical simulations are means of increasing our 
understanding of the physical properties of the undetected 
IGM. Many numerical s imulations of the cosm i c reionisation 
have b een performed ([Nakamoto et alj l200ll; ICiardi et aT 
20031: ISokasian et al.1 120041: llliev et al " 



2006 



et al.l 
20071: 



Trac & Cenl 120071: iFinlator et al .11 20091: lAubert & Tessierl 



20101 : IPetokova fc Springe]|l201ll : lAhnet al.l l2012l ). More- 
over, in order to predict the 21-cm signal emitted in 
intergalactic medium (IGM) during the reionisation epoch, 
which should be observed by radio- interferometers such 
as LOFAR, MWA and SKA, sim ulations of the cosmic 
reionisation have been performed. ([Ciardi &; Madau| [2003: 
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Mell ema et al.l 



Tho mas et al.l 

201 ih . 



20061 : iLidz et a l. 2008: ISantos et all lioOL 
20091 : Baek et al. 2009; Vonlanthe n et al 



The nature of the sources of the cosmic reionisation 
is still uncertain, but stars in galaxies, at least, are ex- 
pected to contribute to the process, since many galax- 
ies have already been observed at z > 6. In a cold 
dark matter (CDM) cosmology, structure formation pro- 
ceeds in a bottom-up manner. Therefore it is expected that 
low mass galaxies were dominant sources of ionising pho- 
tons during the early phase of reionisation. However it is 
also well known that star formation in low mass galaxies 
is very sensitive to feedback effects owing to their shal- 
low gravitational potential. For instance an external ul- 
traviolet background (UVB) radiation field which devel- 
ops as the cosmic reionisation proc e eds heats the gas u 
to - 10 4 K (jThoul fc Weinberg! fl996l: lOkamoto et al.| 12008^ 



and dissociates molecu l es (iRicotti. GnedinTfc Shull 20011 
busa fe Umemural 12004 Hasegawa. Umemura &; Kitavamal 



20091 : I Johnson et al.1 l2012n ~ Moreover ISusal (j2008l ~ have 



shown an external UVB can suppress the star formation even 
in a galaxy that is massive enough to sustain an ionised hot 
gas. The suppression in the massive galaxy is mainly due to 
the fact that photo-heating prevents gravitational inability. 
Not only the external UVB but also UV feedback from its 
own stars is likely to regulate star formation in each galaxy. 
In order to investigate ionising photon esca pe fraction from 
high— z dwarf galaxies, IWise &; Cenl (|2009h have performed 
radiation hydrodynamics (RHD) simulations in which the 
internal UV and SN feedback are included. They have shown 
that the internal feedback evacuates a large amount of gas 
from halos with masses less than 1O 8-5 M0, and regulates 
the star formation rate. If previous episodes of star forma- 
tion and feedback effects associated with the star formation, 
which are not included in their simulations, are taken into 
account, the impact of the internal feedback is expected to 
be more drastic. In addition, such evacuation would lead to 
the enhancement of the escape fraction of ionising photons, 
which is an important factor to regulate reionisation history. 
According to the previous studies, it is strongly expected 
that UV feedback plays a very important role in regulating 
the histories of the cosmic star formation and the reionisa- 
tion. But unfortunately it is often hard to implement the 
internal UV feedback in simulations of the cosmic reionisa- 
tion, since calculation costs for the RT are very high. 

Illiev etHl (|200/t > have explored the reionization history 
regulated by UVB with high resolution iV-body and post- 
processing RT calculations, assuming a constant character- 
istic halo mass of 1O 9 M0 below which stars are forbidden 
to form in ionised regions. However, this characteristic mass 
should depend on the halo formation epo ch, the strength 
of the UVB radiation field, and so on (|Kitavama et al.l 
l200ll : ISusa fe Umemural 12004 : lOkamoto et alJl2008h . More- 
over the internal UV feedback cannot be correctly included 
due to a lack of spatial resolution in RT calculations, al- 
though the spatial and mass resolutions for iV-body sim- 
ulations are sufficiently high to take into account ionis- 
ing photons from small structures. Hence the star for- 
mation history in th eir simulations might be inaccurate. 
iFinlator et al.1 (|201lh have recently performed simulations 
of the cosmic reionisation in which the transfer of UV 
photons is consistently coupled to hydrodynamic calcula- 



tions. Hence they can consider suppression of star forma- 
tion in galaxies due to photo-evaporation without assum- 
ing a characteristic mass. However the spatial resolution 
for their RT calculations is not high enough to correctly 
take into account the internal UV feedback. The lack of 
spatial resolution also implies using an escape fraction pa- 
rameter, which is frequently assumed to be constant over 
all masses of galaxies. Since the ionising photon escape 
fraction likely depends on various properties of the galaxy, 
it is doubtful that the assumption of a constant value is 
valid. In addition, there is little agreement on the value 
of the escape fraction in previous studies (|Gnedin et al 



0; 



l200Sl:IWise fc Cenl2009l :lR azoumov fc Sommer-La esen 2010; 
lYaiima et al.l |201l|) . Hence the employment of the escape 
fraction results in a large uncertai nty of the reionisation 
history. IPetokova &; Springell (|201lh have implemented a 
moment-based RT solver in GADGET, and have simulated 
the cosmic reionisation. With the code, spatial resolutions 
for the RT and hydrodynamics calculations are equivalent 
so that the internal UV feedback can be appropriately con- 
sidered. In addition, they did't have to employ a "galactic" 
escape fraction. However there might be still an unresolved 
interstellar medium (ISM) where some photon absorption is 
expected. Hence they adopted an "ISM" escape fraction to 
take into account the photon losses in the unresolved ISM. 
Although some of their models well produce a plausible his- 
tory of reionisation, the minimum halo mass resolved in their 
simulations, typically 3.0 x 1O 9 /i _1 M but 3.8 x 1O 8 /i _1 M 
in the maximum resolution run, seems too massive to discuss 
the suppression of sta r formation at the epoch of reionisation 
( Oka moto et al.ll2008r ) . Therefore the impacts of the internal 
and external UV feedback on the histories of star formation 
and reionisation have not been well established yet. 

In this paper, we aim to clarify the importances of UV 
feedback on cosmic reionisation history by radiation hy- 
drodynamics simulations. Unlike previous numerical simula- 
tions of the cosmic reionisation, we include non-equilibrium 
chemistry regarding the H2 molecule, which is the most effi- 
cient coolant below 10 4 K in a low-metallicity environment. 
H2 cooling plays a crucial role not only for less massive 
halos with T v i r < 10 4 K but also for massive halos with 
T vir > 10 4 K. Needless to say that halos with T vir < 10 4 K 
cannot collapse without H2 cooling in the early Universe. 
But even in halos with T v i r > 10 4 K, internal structures 
would be affected by thermal processes. We have to trace 
the dynamical and thermal evolution of gas in the cold 
T < 10 4 K regime, since stars are born in cold and dense 
regions. In addition, resolving such internal structures is 
very important to understand how ionisin g photons es- 
cape from a ga l axy (|Fernandez &; ShulfeOllh . Very recently, 
I Johnson etall (|2012n ~ have shown that the global and local 
photo-dissociating radiation regulates cosmic star formation 
history in the epoch of reionization. However the roles of an 
inhomogeneous reionisation process and the internal photo- 
ionisation feedback are still unclear, because they simply as- 
sumed a homogeneous reionisation model. We show in this 
paper not only photo-dissociating radiation but also photo- 
ionisation radiation plays crucial roles on the reionisation 
and star formation history. 

This paper is organized as follow. In Section 2, we de- 
scribe our numerical method and the setup of the simula- 
tions. We show the results of the simulations in Section 3. In 
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Figure 1. Specific ionising photon production rate as a function 
of stellar age. The solid, dotted, and dashed lines respectively 
show the photon production rates of HI, Hel, and Hell ionising 
photons. 



Section 4, we compare the simulation results with observa- 
tions, and also discuss the expected effect of some processes 
that we did not consider in this work. Finally we devote Sec- 
tion 6 to the conclusions. Throughout this paper we assume 
a A cold dark matter universe with cosmological parameters: 
Q b = 0.049, Q = 0.27, Q A = 0.7 3, h 0.71, and g_ 8 = 0.81 
based on 7- year WMAP results (jjarosik et alj|201ll ) 



2 SIMULATION 

In order to self-consistently simulate the reionisation and 
star formation history, we have to couple hydrodynamics to 
the transfer of UV photons from numerous radiation sources. 
For this purpose we use a novel radiation hydrodynamics 
code STA RT (SPH with Tree-based Ac celerated Radiative 
Transfer) ([Hasegawa &; Umemural l2010h , which enables us 
to solve such a numerically expensive problem. We describe 
the detailed simulation method in the following part. 



2.1 Feedback Processes 

2.1.1 Star Formation Recipe 



We adopt the star 



formation recipe employed in 
ISusa &; Umemural (| 20041 ) who explored the impact of UVB 
on the formation of dwarf galaxies. We allow each gas parti- 
cle to form stars, if a gas particle satisfies the following three 
conditions; (i) p > 200p a v, (ii) T < T cr = 5, 000K, and (hi) 
div(v) < 0. The stars are born at a rate of 



dp* 
~dt 



tdy 



(1) 



where p*, p gas , and C* are the star and gas mass den- 
sity, and the star formation efficiency, respectively. The dy- 
namical time tdyn is evaluated as tdyn = 1/ y/kirGp. Here 
the condition (ii) is the most crucial, since this condition 
can be satisfied only if the gas is well shielded against 
UV radiation so that H2 cooling effectively works. How- 
ever the change of T cr hardly affects the final results as 
long as T cr is well smaller than 10,000K, since the cool- 
ing time is well shorter than the dynamical time in a self- 



shielded region (Hasega wa, Umemura &; Kitavamal[2009h . In 
previous simulations of galaxy formation, C* ~ 0.01 is fre- 
quently employed so as to reproduce the Schmidt-Kenn icutt 
relation (|Navvaro fe Steinmetd I2OO0I ; ISaitoh et"alll2008h . In 
addition, Susa & Umemura (2004) have shown the final stel- 
lar mass fraction in a halo is almost independent on the value 
of C* in the case where the star formation is regulated by 
UV feedback, although the star formation history is slightly 
delayed for smaller value of C*. Therefore we carry out all 
of the simulations with C* = 0.03. 



2.1.2 SN feedback 

We include thermal feedback by SN explosions. Firstly, in 
advance of the s imulations, we compute populati on synthesis 
using PEGASE (jFioc & Rocca-Volmerange|[l997h . Assuming 
a Salpeter IMF with masses of IMq-IOOMq and an instan- 
taneous starburst, we obtain a specific SN rate as a function 
of a stellar particle age. According to the rate, we feed the 
neighboring 50 particles around the star particle with SN en- 
ergy in thermal form, assuming a released energy of 10 51 erg 
for one SN explosion. 

Chemical feedback is not included in this study. The in- 
clusion of metal enrichment enhances cooling rate. In addi- 
tion, dust grains likely play an important role; H2 formation 
on the surface of dust grains would increases H2 fraction 
and the absorption of UV photons by the grains leads to 
a decreased UV photons escape fraction. We will include 
chemical feedback in a forthcoming paper. 



2.1.3 UV Feedback 

We use the RHD code START. Here we briefly describe 
the RT method implemented in the code. The implemen- 
tation of ray-tracing in START is basically the same as that 
adopted in RSPH, which has been used for simulations of 
struc t ure formation in the early Universe (ISusa &; Umemural 
2006; Susa 2007, 2008; Hasega wa, Umemura fe Susall2009l : 
ISusa, Umemura &; Hasegawa! 120091 ) . Since the rays are cast 
from all radiation sources to all SPH particles, the RSPH 
scheme allows us to solve RHD problems very accurately 
(jlliev et al.|[2006al : ISusal 120061 : llliev et alJ l2009h . Moreover, 
this ray-tracing algorithm enables us to solve RT without 
losing the spatial resolution of hydrodynamics calculation, 
since SPH particles are directly used for RT grids. With 
RSPH, however, the calculation cost is proportional to a 
number of radiation sources, and it becomes difficult to solve 
problems in which numerous radiation sources should be 
considered. In order to avoid this difficulty, a tree-based ac- 
celeration algorithm is employed in START. In the tree-based 
acceleration algorit hm, we first make an oct-tree structure 
(Bar nes &; Hut|[l986l ) for radiation sources. To make the oct- 
tree structure, we create the root cell which contains all 
radiation sources, and recursively divide the cell into eight 
cells down to the level where each cell contains one radia- 
tion source. Then we check if a cell is distant enough from 
a target SPH particle by using following criteria; 



(2) 



where d s is the distance between the cell and the target 
SPH particle, l s is the size of the cell, and # C rit is the tol- 
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Table 1. List of runs 



Name 



box ; 



[Mpc 



Numbers of particles RT photodissociation f e ism 



rN256L5 


5 


256 3 


X 


2 


no 








rN256L5RT 


5 


256 3 


X 


2 


yes 


shielding 


5 function 


1.0 


rN128L5 


5 


128 3 


X 


2 


no 








rN128L5RT 


5 


128 3 


X 


2 


yes 


shielding 


5 function 


1.0 


rN128L5RT-fe08 


5 


128 3 


X 


2 


yes 


shielding 


5 function 


0.8 


rN128L5RT-fe06 


5 


128 3 


X 


2 


yes 


shielding 


5 function 


0.6 


rN128L5RT-fe04 


5 


128 3 


X 


2 


yes 


shielding 


5 function 


0.4 


rN128L5RT-fe02 


5 


128 3 


X 


2 


yes 


shielding 


5 function 


0.2 


rN128L5RT-noLW 


5 


128 3 


X 


2 


yes 


no 


1.0 


rN128L5RT-nfsh 


5 


128 3 


X 


2 


yes 


optically thin 


1.0 



z=1 1 .1 



s o 



s o 



s o 



-2 





X [Mpc] 



z=7.58 




z=5.99 



7 

1 ? 



1 



1 



X[Mpc] 



X [Mpc] 



Figure 2. Maps of physical quantities in a slice cut through the mid-plane of the simulation box at redshifts z = 11.1, 7.58, and 5.99 
from left to right. The evolution of comoving number density, temperature, and ionisation degree is shown in the upper, middle, and 
lower rows, respectively. 



erance parameter. If the above condition is satisfied, the 
radiation sources in the cell are regarded as one luminous 
virtual source. Adopting this algorithm, the calculation cost 
is dramatically reduced. We have to adopt a sufficiently 
small value of # C rit to obtain accurate results. fl C rit = 0.5 
is a reasonable cho ice (Hasega wa &; Umemural ([2010) , see 
also lOkamoto etafl (|201lh ). 



We consider the photo- ionisation of HI, Hel, and Hell 
by solving the transfer of ionising photons appropriately. 
The cross-sections are respectively given by 



cr m {v) = 6.30 x 10~ 18 (^J 



(3) 
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cthiO) = 7.42 x 10" 



crHeiiO) = 1.575 x 10" 18 



cm , 
(4) 

(5) 



where z/i, and are the Lyman limit frequencies of 
HI, Hel, and Hell, respectively. The corresponding photo- 
ionisation rates are expressed as 



oexp(-r^) 



hv 



o exp(-T„) 
hv 



<JYL\{y)d£ldv 



aHel(^)rfQdz/, 



-CTHeIl(V)G^G^. 



(6) 



(7) 



(8) 



As shown in previous studies (iTenorio-Tagle et"aT] 1986; 
iFrenk Mellemal 1 19941 : iNakamoto et al.l l200lh . integrated 
photoionisation rates can be obtained knowing the optical 
depths at the three Lyman limit frequencies only, since we 
already know the shape of intrinsic intensity i^o and the 
dependence of the cross-sections on frequency. We produce 
a three dimensional, i.e., a function of r Ul , 7v 2 , and r U3 , 
look-up table by integrating photo-ionisation rates over all 
frequency ranges as a preprocessing of the simulations. By 
using the look-up table we can appropriately take into ac- 
count the dependence of their cross-sections on frequency. A 
similar method is applied to calculate photo-heating rates. 

Hydrogen molecules are easily dissociated by Lyman- 
Werner (LW) band (11.2-13.6eV) photons via the two 
step Solomon process. We take the photo-dissociation of 
H2 into account by calculating the column density of 
H2 and adopting the s elf-shielding function introduced by 
iDraine fe Bertoldil {l996). Then the flux in the LW band 
can be written as 



Flw — Fuw,ofs (A^H 2 ,14) 



(9) 



where Flw,o is the incident flux, A^h 2 ,i4 is the H2 column 
density in units of 10 14 cm -2 , and 



1, 



"3/4 



X < 1 
X > 1 



(10) 



We also include photo-detachment of H~ and photo- 
dissociation of H 2 , assuming these species are optically thin 
against the radiation owing to their considerably smaller 
abundances. These rates can be expressed as 



L J hv 



dQdv, 



(ii) 



where hv H - = 0.74eV and hv H + — 0.062eV. The cross- 
sections are taken from iTegmark et al.l (|l997h and IStancill 
(|l994l ). respectively. Here we should note that the self- 
shielding function for H2 photo-dissociation was derived 
with the assumption that a medium is static. Consequently 
the shielding effect might be stronger than in the case where 
velocity gradients are properly included. Therefore we addi- 
tionally simulate two extreme cases, by changing the treat- 
ment of H2 photo-dissociation; (i) a run in which we never 
use the shielding function ( H2 molecules are optically thin 



in the LW band), (ii) a run in which photo-dissociation pro- 
cesses (not only H2 photo-dissociation but also H~ photo- 
detachment and H^ photo-dissociation processes) are ex- 
cluded. We verify the importances of the photo-dissociation 
processes by performing these two runs. 

We compute the HI, H el and Hell ionising photon pro - 
duction rates in PEGASE (jFioc Rocca- Volmerangd[l997n , 
in the same way that we evaluate the specific SN rate as a 
function of stellar age. We plot the photon production rates 
as a function of stellar age in Fig[T] Here, for simplicity, 
we assume that the shape of spectral energy distribution 
(SED) always corresponds to that of a black body spectrum 
with an effective temperature of 50,000K. The assumption 
is reasonable due to the following reasons; (i) UV radiation 
is dominantly emitted by young stars, (ii) With the IMF 
we choose, the SED of young stellar population can be well 
approximated b y a black body w ith an effective tempera- 
ture of 50,000K (jBaek et alJ l20Q9). Evaluating reaction rates 
between UV radiation and matters by the way described 
above, we implicitly solve a non-equilibrium chemistry re- 
garding nine species, i.e., e~, H + , H, H~, H2, Hj~, He, He + , 
and H e 2+ . The initial fractions are taken from Gal li &; Pallal 
(| 19981 ). 

Finally, there is a synergy effect between SN and UV 
feedback. It is often true that thermal feedback by SN explo- 
sions cannot work efficiently, because radiative cooling pro- 
cesses efficiently work in a high density regions, which are 
the sites of star formation, and gas in these regions turns out 
to be cooled quickly. In our simulations, a timestep is basi- 
cally much smaller than the lifetime of massive stars. Hence 
UV feedback works prior to the first SN explosion. The UV 
radiation from newly born stars heats the adjacent medium, 
and the gas density decreases before thermal energy released 
by SN explosions is poured into the gas. Consequently, SN 
thermal feedback works more efficiently when UV feedback 
is included. 



2.2 Setup 

Each simulation starts f rom z — 150 with i nitial conditions 
generated with Grafic-2 ([Bertschingerl I200H ) . We performed 
simulations with two different resolutions. In high resolution 
runs, we handled 256 3 particles each for SPH and DM parti- 
cles in a volume of comoving Lbox = 5Mpc on a side. In this 
situation, the particle masses correspond to 4.9 x 1O 4 M0 
for SPH and 2.5 x 1O 5 M for DM, respectively. Hence a 
small galaxy of Mh a io,min ~ 2.5 x 1O 7 M can be resolved 
with 100 particles. We note that this mass resolution is 
~ 20 times better than that in the highest resolution run 
of IPetokova &; Springell (|201lh . The gravitational softening 
length is set to 0.03 x Lbox/Af 1//3 , which corresponds to 595 
comoving pc. In low resolution runs, on the other hand, we 
handled 128 3 particles each for SPH and DM particles in 
the same size box. Hence the minimum halo mass in the low 
resolution run is Mh a io,min ~ 2.0 x 10 8 Mq which is simply 
8 times higher than that in the high resolution run. Since 
the spatial resolutions of both the RT and hydrodynamics 
calculations are sufficiently smaller than the typical size of 
small galaxies, we can appropriately take into account in- 
ternal UV feedback in each galaxy. The spatial resolution of 
the high resolution RT simulation is roughly equivalent to 
those of previous numerical works that have evaluated ionis- 
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ing photo n escape fractions of ionising photons from galac- 
tic h alos (|Gnedin et al.l [20081 ; iRazoumov &; Sommer-Laesenl 
2010). Although photon losses in galaxies can automatically 
be taken into account owing to the high spatial resolution, 
there still might be absorp t ions in the spatially unresolved 
ISM. IPetokova fe Springell (|201lh have adopted an "ISM" 
escape fraction / e ,isM to account for such unresolved absorp- 
tions. We note that this "ISM" escape fraction has a com- 
pletely different physical meaning from the "galactic" escape 
fraction which has be en usually utiliz e d in previous simula - 
tions on reionisation (|lliev et alJl2007l : iFinlator et al.ll201ll ). 
In our simulations we basically work with / e ,isM = 1.0, but 
we also carry out simulations with lower values of the ISM 
escape fraction. We should emphasize that the variance of 
galactic escape fractions connected to the different proper- 
ties of each halo is automatically taken into account even 
if we employ an ISM escape fraction. The validity and de- 
pendence of our results on the ISM escape fraction will be 
discussed in Section [4721 Parameter sets for all runs are sum- 
marized in Tab lei. 



3 RESULTS 

3.1 Global Evolution 

In this section, we show the evolution of global quantities, 
such as ionised fraction and star formation rate density. We 
would like to point out that our simulation box is too small 
to mitigate the effect of cosmic variance on these quantities. 
Although the presented global quantities cannot be consid- 
ered to be representative of the whole Universe, these quan- 
tities are very useful to understand the roles of UV radiation 
feedback. In addition, these quantities can be used for the 
comparison with the results of similar studies. To visually 
understand a simulated reionization process, we show the 
simulated maps of hydrogen number densities, gas temper- 
atures, and ionization degrees at redshifts of 11.1, 7.58, and 
5.99 in Fig. [2] This is for the rN256L5RT run, but a similar 
evolution can be seen in other runs. After stars are born at 
the early epoch of z ~ 11, well known "butterfly-shaped" 
ionised regions associated with star forming regions appear. 
Each ionized region grows, and the half of the simulation 
volume is ionized at z ~ 7.5. The reionization is finally com- 
pleted by z ~ 6. 

In the top panel of Fig(3] we show the simulated ionisa- 
tion histories. As we can see from the figure, in both cases 
of the high and low resolution run, the mass weighted aver- 
age ionisation fractions of hydrogen are slightly higher than 
volume weighted ones. This trend is indicating that high 
density regions are reio nised prior to lower density regions , 
as in previous studi e s (llliev et a . 1 l2007l : iTrac fc CenN2007l : 
IFinlator et all 120091 : iBaek et al.l l2009h . Obviously, ionised 
fractions at high redshifts are higher in the high resolution 
run, since the first UV source tends to appear at higher red- 
shift as the mass resolution gets higher. 

In the bottom panel of Fig. [3] we represent the simu- 
lated star formation rate densities (SFRDs). We also plot 
the SFRDs inferred from observed lumi nosity functions of 
high-;? drop-out gal axies in the figure (lOuchi et al.l [2009; 
iBouwens et al.ll201lh . It has often been argued that stars are 
overproduced in numerical simulations. Actually our simu- 
lations also seem to overproduce stars, if we exclude UV 



rN256L5RT <mass weighted> 
rN256L5RT <volume weighted: 
rN128L5RT <mass weighted> 
rN128L5RT <volume weighted: 




Figure 3. Upper panel: Hydrogen ionisation histories. The solid 
and dashed curves respectively indicate the mass and volume 
weighted ionisation histories in the rN256L5RT run. while the 
dotted and dat-dashed curves respectively indicate the mass 
and volume weighted ionisation histories in the rN128L5RT run. 
Lower panel: Cosmic star formation rate densities. The cos- 
mic star formation rate densities in the rN256L5RT, rN256L5, 
rN128L5RT and rN128L5 runs are respectively shown by the thick 
solid, thick dashed, thin solid, and thin dashed curves. The filled 
and open squares are respectively the extinction corrected SFRD 
integrated down to L = 0, and the no extinction corrected SFRD 
integrated down to L « 0.1L* from Ouchi et al. (2009). The filled 
circles are the SFRD integrated down to 0.05L* =3 from Bouwens 
et al. 2011. 



feedback. It is interesting that the SFRDs in both runs are 
suppressed by a factor of ~ 5 — 6 at z = 6 if we include 
UV feedback. This strong suppression makes the simulated 
SFRDs match the SFRDs inferred by observations at z ^ 8. 
At z rsj 10, the simulated SFRDs remain inconsistent with 
the observation. However whether simulations and obser- 
vations are consistent or not might be inconsequential at 
this stage, because there are some uncertainty in the both 
of the numerical and observational studies. Concerning our 
study, as mentioned above, the volume size of the simula- 
tion is too small to be regarded as a representative of the 
Universe. Moreover some processes expected to be impor- 
tant for the star formation process, that will be discussed in 
a later section, haven't been included yet. As for observa- 
tional studies, on the other hand, the shape and amplitude 
of the luminosity function of faint galaxies are still uncer- 
tain especially at high- z . In addition, the fact that our small 
simulation box does not contain luminous galaxies that are 
directly observed illustrates the mismatch between our sim- 
ulations and the observations. Thus here the point to which 
we should pay attention is that UV feedback significantly 
affects the cosmic star formation history. C omparing the 
results of similar previous studies, such as Finl ator et al.l 



The Impacts of UV feedback on galaxies during the EoR 7 



■ '\ 

\ 

\ 

\ \ 


J 




/ \ rN256L5RT 

'i^l256L5 

rNTi28L5RT 

^ rN1^5 




\ 





















Mhalo/Msun 

Figure 4. Baryon mass fractions (upper panel) and star forma- 
tion rates (lower panel) at z = 6.0 as a function of halo mass. The 
curves are mean values at a given mass range, and the shaded re- 
gions attached to the curves are their dispersions. The thick solid, 
thick dashed, thin solid and thin dashed curves respectively show 
the results by rN256L5RT, rN256L5, rN128L5RT and rN128L5. 
The thick and thin dot-dashed lines in the upper panel indicate 
the critical baryon mass fraction given by Eq. 1121 Halos located 
above this lines contain more than 100 baryon particles. 



(|201lh and IPetokova fe Springell (|201lh . we find that the 
suppression seen in our simulation is considerably stronger 
than those shown in these studies. We should mention that 
the former have not appropriately considered the internal 
UV feedback, while the latter do not have a sufficient mass 
resolution to discuss the suppression. Also none of these two 
studies take into account H2 chemistry. We study the reasons 
for our strong suppression of star formation in the following 
sections. 

In the lower panel of Fig. [3l we can see that the SFRD 
in rN256L5RT at z < 9 is lower than that in rN128L5RT, 
in spite of the fact that SFRDs in the corresponding runs 
without UV feedback, namely rN256L5 and rN128L5, are 
almost identical at this epoch. Relating to the discrepancy 
of the SFRDs, the ionization history at the later stage in 
rN256L5RT is delayed compared to that in rN128L5RT, as 
shown in the upper panel of Fig. [3] These discrepancies are 
thought to be caused by numerical resolution effects. In Sec- 
tion [43] we discuss it in detail. 

3.2 Properties of Halos 

In the previous subsection, we showed that the cosmic star 
formation rate is strongly suppressed by UV feedback effects. 
In this section we analyze some properties of the halos to 
uncover the reasons for the suppressed star formation. To 
identify individual halos, we first use the friends-of-friends 



method with a linking length of 0.2 x N 1 ^ 3 . Then, gradually 
shrinking a sphere of given radius, we refine the location 
of the centre of a halo to be the mass centre of it. When 
the density within the sphere is equal to the virial density, 
we define the radius as the virial radius of the halo. We 
define the ratio of the amount of baryon particles within a 
halo's virial radius to the total halo mass as the baryon mass 
fraction of the halo. 

Since the suppression of the star formation is the most 
noticeable at z ~ 6.0, we here focus on halos at this epoch. In 
the upper panel of Fig|4] we compare the baryon fractions in 
the rN256L5RT, rN256L5, rN128L5RT, and rN128L5 runs. 
Before we show the results in detail, we should mention 
the spatial resolution of the presented simulations. As de- 
scribed above, before the feedback processes start to work, 
there is a sufficient number of SPH particles to resolve the 
internal feedback in each halo. However, once the internal 
feedback starts to work, SPH particles are evacuated from 
halos. As a result, the spatial resolutions of the hydrody- 
namics and RT calculations, especially for low mass halos, 
gradually decreases as time goes by. As shown by previous 
studies, at least 100 SPH particles are required to resolv e 
each halo (jPetokova fe Springell lioill ; iFinlator et aUfeOllh . 
The baryon mass fraction which satisfies this critical condi- 
tion can be given by 



fb,c 



5.4 x 10" 



10 8 M< 



(12) 



for the high resolution runs. The corresponding equation for 
the low resolution runs is easily obtained by multiplying the 
above equation by eight. We trace the corresponding lines in 
the upper panel of FigU Then we can see that the internal 
structures of halos with M ha io < 2 x 1O 8 M (or M ha io < 
7 x 10 8 M© in the low resolution run) cannot be discussed 
with our simulations. Indeed the baryon mass fractions of 
halos with M ha io = 2 - 5 x 10 8 in rN128L5RT show larger 
deviations than those in rN256L5RT does. It is thought to 
result from a very small number of gas particles in the halos. 
Therefore we should consider the internal quantities of halos 
less massive than these mass scales to be unreliable. 

Roughly speaking, the baryon mass fractions of halos 
with Mhaio > 1O 9 M0 are insensitive to UV feedback, though 
the baryon mass fractions are slightly smaller than the value 
expected from the cosmological parameters Q D /^o — 0.18. 
On the other hand, the discrepancy of baryon mass fractions 
between the runs with and without UV feedback is remark- 
able for halos with M ha io < 1O 9 M . At z = 6.0, the char- 
acteristic mass scale below which photo-evaporation effec- 
tivel y takes place is roughly consistent with the value shown 
bv lOkamotoeTall diooi ). although they assume an homo- 
geneous external UV feedback and do not consider inter- 
nal UV feedback and self-shielding. This consistency likely 
comes from the following facts: (i) the temperature of ionised 
gas is not changed dramatically even if internal UV feedback 
is included, (ii) for less massive halos, self-shielding is no t 
as efficient as for massive halos ([Taiiri &; Umemu ra 1998). 
Therefore, as far as we consider the epoch when the Uni- 
verse is almost ionised, the characteristic mass investigated 
with a uniform UV background intensity would be approxi- 
mately correct. But we should mention here that the internal 
structure of halos is significantly affected by the internal UV 
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Figure 5. Maps of the distributions of H nucleus number density (right panels), neutral H number density (central panels), and H2 
fraction(left panels). The upper row is for Halo-R (a massive halo in rN256L5RT run), while the lower row is for Halo-H (the same halo 
but without UV feedback). 



feedback even in massive halos as we will show in Section 

E23I 

We represent the star formation rate as function of 
halo mass in the lower panel of Fig[4] As expected from 
the top panel of Fig[4j the star formation rates of halos 
with Mhaio < 1O 9 M is strongly suppressed via photo- 
evaporation. An interesting result can be seen in the star 
formation rates of massive halos. The star formation in such 
massive halos is also suppressed, despite the fact that the 
massive halos hardly lose any mass. In the following two 
subsections, we show reasons for this phenomenon. 



3.2.1 Impact of UV feedback on Internal Density 
Structures 

A possible mechanism which causes the suppression is the 
decrease of gas density driven by the internal feedback. Soon 
after the birth of stars in a dense cooled region, the stars 
emit ionising photons which quickly heat up the surrounding 
gas up to rsj 10 4 K. This process leads to the destruction of 
the dense star- forming region by the high thermal pressure. 
Also the propagation of the D-type i onisation front makes 
blow-out by SN explosions effective (IKitavama Yoshidal 
2005). According to the equation (1), the star formation 
rate is proportional to p 1 ' 5 . Therefore if the gas density is 
decreased, this directly indicates that the further star for- 
mation in the region is suppressed. 

We choose one massive halo with M to t ~ 6.7 x 1O 9 M 
from rN256L5RT. We also pick the same halo from rN256L5, 
and compare them in Fig. Hereafter we call the former 
Halo-R, and the latter Halo-H. In each halo, more than 



20,000 SPH particles are included. Thus we can discuss 
their internal density structures. As we can see from Fig. \5\ 
the internal density structures of these two halos are quite 
different. Some clumpy structures can be seen in Halo-H, 
that are significantly smoothened by the UV feedback in 
Halo-R. To emphasize the difference in the internal den- 
sity structures more quantitatively, we show cumulative his- 
tograms of the mass fraction where the overdensity S ex- 
ceeds a given value in Fig[6j Here we define the overdensity 
as 5 = p h / < p h >haio, where < p h > ha io is the average gas 
mass density of the massive halo. As shown in Fig.[6j Halo-H 
has more high density components than Halo-R. This is the 
direct cause for the suppressed star formation in the massive 
halo. The destruction of clumpy structures likely works more 
effectively in less massive halos. Since a massive halo forms 
through coalescences of smaller halos which are sensitive to 
the UV feedback in a CDM cosmology, resolving previous 
episodes of star formation would leads to stronger suppres- 
sion of the star formation. Indeed, in our simulations, the 
sensitivities of the star formation in massive halos to the 
UV fe edback are higher than those shown by IWise &; Cenl 
(l2009h. 



It is also expected that the internal UV feedback af- 
fects the ionizing photon escape fraction from a galaxy, since 
the escape fractions strongly depend on the internal struc- 
ture of a gala xy, such as clumsiness and the covering factor 
of clumps ([Fernandez &; Shullll201lh . In addition, stars are 
usually born in a high density region where a high recombi- 
nation rate is expected. Hence if we don't take into account 
the back-reaction by UV radiation emitted by newly born 
stars, the ionising photons hardly escape from the star form- 
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Figure 6. Fractions of baryon mass that resides in a place where 
the baryon overdensity 5 exceeds a given value. The solid and 
dashed curves represent Halo-R and Halo-H respectively. 



ing regions. In lUmemura et ah 1 (|2012T ). we have explored the 
difference between the ionizing photon escape fractions from 
Halo-R and from Halo-H, by solving the RT of ionising pho- 
tons in post-processing. As a result, we have found that the 
escape fractions of the HI, Hel, and Hell ionising photons 
from Halo-R are respectively 0.31, 0.35, and 0.13. These val- 
ues were indeed higher that those from Halo-H by a factor 
of 2. 

Our mass resolution does not allow as to resolve the 
structure of the low mass clouds where stars are born. There- 
fore we should consider the possibility that the results pre- 
sented in this section would be altered if we increase the nu- 
meri cal resolution. R ecent high resolution RHD simulations 
by I Dale et al.l (J2012|) have shown that the dynamical effect of 
photo-ionisation disrupts low-mass star-forming clouds with 
1O 4_5 M0 within a few Myr. This result indicates that fur- 
ther star formation in these regions is suppressed and the 
feedback enhances the escape fractions of ionizing photons 
from the star- forming clouds. Thus we expected that the 
results in this section are qualitatively unaffected by the 
numerical resolution, even though the quantitative values 
might be affected. 



3.2.2 Dissociation 0/H2 molecules 

In previous studies of reionisation by numerical simulations, 
the non-equilibrium chemistry of H2 has often been ne- 
glected for simplicity. As mentioned above, if we want to 
resolve the internal structure of halos, which is closely re- 
lated with the star formation rate and with the escape prob- 
ability of ionising photons, we have to consider cooling pro- 
cesses below 10 4 K. In the low metallicity environment which 
is prevalent in the early Universe, H2 molecules are the most 
efficient coolants in the low temperature regime. Therefore, 
if the H2 fraction is considerably decreased, the gas cannot 
cool down to T < 10 4 K and stars are hardly formed. 

In this section, we show how severely the H2 fraction 
of each halo is affected by UV feedback, i.e., the photo- 
ionisation and photo-dissociation processes. Before we show 
the results in detail, we briefly introduce some processes 
closely related to the determination of H2 fraction before 
we discuss the results. The destruction of H2 molecules is 
caused not only by photo-dissociating radiation but also 
indirectly by photo-ionising radiation. In the photo-ionised 
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Figure 7. Mass weighted average H2 fractions in high density 
regions (5 ^ 1) of halos at z = 6 as a function of halo mass. 
The dotted, dashed and sold lines show the results of rN128L5, 
rN128L5RT-noLW and rN128L5RT, respectively. The vertical 
dot-dashed line indicates the mass of halo which typically includes 
100 SPH particles in the rN128L5RT run. 
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Figure 8. Ionised hydrogen fractions (upper panel) and tem- 
peratures (lower panel) in high density regions (5 ^ 1) of halos 
at z = 6 as a function of halo mass. The dotted, dashed and 
sold lines show the results of rN128L5, rN128L5RT-noLW and 
rN128L5RT, respectively. The vertical dot-dashed line indicates 
the mass of halo which typically includes 100 SPH particles in 
the rN128L5RT run. 



hot regions, H2 molecules are collisionally dissociated. More- 
over, the smoothing of internal density structures shown in 
the previous section weakens the self-shielding effect. On the 
other hand, there is a process which promotes H2 formation. 
In the absence of dust grains, H2 molecules are produced via 
the H _ process in which electrons work as catalysts. Thus 
as the e~ fraction increases, H2 formation rate also increases 
with a given temperature. Actually, the promotion of H2 for- 
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mation always works when an enhancement of the electron 
fraction exists, but usually the positive feedback is remark- 
able only if the photo-ionisation rate is not too high and 
the medium is only slightly ionized. Otherwise the strong 
negative photo-ionisation feedback becomes dominant. 

As already mentioned, the internal quantities of low 
mass halos are not reliable due to the very small number of 
gas particles in the halos. Therefore we will only discuss the 
internal quantities in massive halos with masses ^ 1O 9 M0 
about which we are currently curious. In addition, we here 
focus only on high density regions, since stars are expected 
to form there. We show the mass weighted average H2 frac- 
tion in the high density regions of halos at z = 6.0 in Figl?] 
To emphasize the effects of photo-ionisation, we also plot 
the mass averaged ionised fraction and temperature in high 
density regions in Fig. [8] We define high density regions as 
the places where 5 ^ 1.0, where the definition of 5 is the 
same as in the previous section. We show not only the H2 
fraction in the rN128L5 and rN128L5RT runs but also in 
the rN128L5RT-noLW run, in which the photo-dissociation 
and photo-detachment processes are excluded, in order to 
evaluate impact of the photo-dissociating radiations. 

In the massive halos, the H2 fractions in both the 
rN128L5RT and rN128L5RT-noLW runs are considerably 
smaller than in the rN128L5 run. Moreover, the differ- 
ence in the H2 fractions between in the rN128L5RT and 
rN128L5RT-noLW runs is not so large. This result tells us 
two facts. First, star formation in massive halos is inhibited 
by the destruction of H2 molecules. Second, the destruction 
of H2 is caused by photo-ionisation(-heating) more than by 
photo-dissociation. As obviously shown in the upper panel 
Fig. the mass weighted ionised fraction tends to increase 
as the halo mass increases, though the ionisation degrees 
are rather small. This trend basically originates in the fact 
that the specific star formation rate becomes higher as the 
halo mass increases, as shown by the bottom panel of Fig 
ID A similar trend can also be seen in the average tem- 
peratures of halos. The gas in massive halos is heated up 
to - 8,000 - 10,000K, in spite of the fact that they are 
not so highly ionised. It has been known that such a mild 
photo-heating is caused by high energy photons. High en- 
ergy photons can travel longer paths than lower energy pho- 
tons, since the cross-sections for photoionisation processes 
become smaller as the energy of incident photons increases 
(equations 3, 4, and 5). In addition, the heating rate by one 
high energy photon is higher than that by a lower energy 
photon, since the released electron has a larger kinetic en- 
ergy. As a result, the gas is mildly heated up even if it is 
almost neutral. Since star formation in massive halos is con- 
tinuous, the average H2 fractions in these halos would take 
an almost-equilibrium value. Indeed, the average H2 frac- 
tions in the massive halos, about 10 ~ 5 — 10 ~ 4 as shown in 
Fig. [7] are roughly consistent with the equilibrium value at 
- 8,000 - 10,000K with / H n « 0.1 - 0.2. Here we should 
mention that such mildly heated regions are thought to be 
the places where the star formation is originally possible in 
rN128L5 run. It is also true that photo-dissociating radiation 
can affect higher density regions than photo-ionising radia- 
tion can, but the amount of gas residing in the regions solely 
affected by the photo-dissociating radiation are smaller than 
that in the mildly heated up regions. Hence we conclude that 
the suppression of star formation in massive halos is mainly 
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Figure 9. Star formation rates of halos at z = 6 shown as a 
function of halo mass. The solid, dashed, dotted, and dot-dashed 
curves indicate the results in the rN128L5RT, rN128L5, rN128L5- 
noLW, and rN128L5-nfsh runs, respectively. The vertical dot- 
dashed line indicates the halo mass which typically includes 100 
SPH particles in the rN128L5RT run. 

caused by the mild photo-heating and the resultant colli- 
sional dissociation in addition to the destruction of clumpy 
structures shown in the previous section. 

If we focus on halos with masses ~ 5 x 1O 8 M0, we notice 
an interesting fact that the H2 fractions in the rN128L5RT- 
noLW run are slightly higher than those in the rN128L5 run. 
This trend might result from the positive feedback due to the 
enhancement of the electron fraction, but we consider these 
halos not to be sufficiently resolved to draw this conclusion 
with any certainty. 



4 DISCUSSION 

4.1 Dependence of the Results on Treatment of 
Photo-dissociation 

As we already mentioned in Section [2.1.31 the self-shielding 
function we employ might lead to a stronger shielding ef- 
fect than if we correctly solve the transfer of LW band pho- 
tons. Hence we should analyse how our results change with 
the treatment of the photo-dissociating process. In Figl9j we 
show the star formation rate of halos. As simply expected 
from Fig. in massive halos, the star formation rates in 
rN128L5RT and rN128L5RT-noLW are nearly equivalent. 
On the other hand, for less massive halos, the star forma- 
tion rates in rN128L5RT-noLW are slightly higher than that 
in rN128L5RT due to the larger fraction of H2 molecules. 
We also plot the star formation rates of the rN128L5RT- 
nfsh run. Since we make the optically thin approximation 
for LW radiations in the rN128L5RT-nfsh run, the photo- 
dissociation process is maximal in this run. We find that the 
discrepancy between the star formation rates in rN128L5RT 
and rN128L5RT-nfsh runs is a factor of ~ 1.5 — 2.0 at most. 
We also show the resultant reionisation histories in Fig llQl 
We find that the difference of the overlapping redshifts, when 
the average neutral fraction is 0.5, in the rN128L5RT and 
rN128L5RT-nfsh runs is just Az « 0.4 corresponding to 
Ar e ~ 0.036. Therefore we expect that our results would 
not change dramatically if we solved correctly the transfer 
of LW photons. 
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Figure 10. Evolution of volume weighted average HII fractions 
(top panel) and stellar mass fractions (lower panel) shown as a 
function of redshift. The solid, dashed, dotted, and dot-dashed 
curves indicate the results in the rN128L5RT, rN128L5, rN128L5- 
noLW, and rN128L5-nfsh runs, respectively. 



4.2 HI Fraction at the End of Reionzation and Its 
Dependence on the Choice of ISM Escape 
Fraction 

In this study, we never assumed a "galactic" escape frac- 
tion, since our calculation automatically takes the absorp- 
tion within a halo into account. We also assumed no photon 
losses in the subgrid, unresolved medium, and so worked 
with an ISM escape fraction of / e ,isM = 1.0. But unfortu- 
nately there is no guarantee that our simulations sufficiently 
resolve absorption in the ISM. Hence we first show how our 
results change with different choices of / e ,isM. We show ion- 
isation histories simulated with / e ,isM = 0.2 — 1.0 in the 
upper panel of Fig llll As shown in the figure, ionisation 
history strongly depends on the choice of / e ,iSM- Obviously, 
reionisation proceeds more quickly as / e ,isM increases. Here 
we should point out an interesting fact; in spite of the strong 
dependence of the reionisation history on / e ,isM, the star for- 
mation history is almost independent on / e ,isM as we can see 
in the bottom panel of Fig. 1111 This implies that the global 
star formation at such high redshifts is mainly regulated 
by the internal feedback processes rather than the external 
UVB, at least in our simulations. Note however that it might 
not be always true, since our simulation volume i s too small 
to be r epresentative of the Universe. As shown in Iliev et al. 
(2006b), the global reionisation and star formation histories 
vary considerably such small volumes. Moreover we cannot 
consider UV radiation emitted from distant sources properly 
because of the limited volume. 

In Fig |12l we plot the simulated neutral hydrogen frac- 
tions and compare them with neutral hydrogen fractions ob- 
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Figure 11. Evolution of the volume weighted average HII frac- 
tion (upper) and stellar mass fraction (lower). The results com- 
puted with / e ,isM — 1-0 (solid curve), / e ,iSM — 0-8 (dashed 
curve), / e ,isM — 0-6 (dotted curve), / e ,lSM — 0-4 (dot-dashed 
curve), and / e ,iSM = 0.2 (dot-dot-dashed curve) are shown. 



servationally estim ated from the spectra of 19 QSOs at z ~ 6 
(|Fan et al.ll2006bh , from Ly q damping wing of Gamma-ray 
burst 050904 at z = 6.3 (jTotani et all l2006h . and from 
the evolution of t he LAE luminosity function at z = 6.6 
([Ouchi et aLl feoiO). As shown in the figure, our simulation 
results and observed neutral hydrogen fractions are well con- 
sistent if we assume / e ,isM greater than ^0.5, otherwise sim- 
ulated neutral fractions are larger than the observational val- 
ues. The ISM escape fraction should approach unity as the 
mass resolution increases. Unfortunately it is still unclear 
how much resolution we need to remove such an uncertainty. 
But we think / e ,isM > 0.5 for our s imulations is reason able, 
because of the following arguments. I Wise &: Cenl (J2009) have 
shown that an ionizing photon escape fraction from a galaxy 
with mass of ~ 1O 8 M0, which roughly corresponds to the 
minimum halo mass in our low resolution run (rN128L5 se- 
ries), is 0.4 — 0.7. The maximum spatial resolution in their 
simulations is O.lpc so that the dumpine s s in a galaxy is 
possibly sufficiently resolved. IWise &; Cenl (|2QQ9h have also 
commented that the computed escape fractions should be 
considered as lower limits, since their simulations have not 
taken into account previous episodes of star formation which 
likely decreases the baryon mass fraction as shown in this 
paper. In addition, needless to say that the galactic escape 
fraction should be smaller than the ISM es cape fraction since 
the former includes the latter. Moreover, Dale et al.l (|2Q12h 
have shown the escape fractions from clouds with 1O 4_5 M0, 
which roughly corresponds to the mass resolution of our sim- 
ulation s, are typically gr e ater t han 0.6. 

In lUmemura et al. I (|2Q12h , we have computed the es- 
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Figure 12. The evolution of the volume and mass weighted aver- 
age HI fractions is shown in the upper and bottom panels, respec- 
tively. The results computed with / e ,lSM = 1-0 (solid curve), with 
Je,lSM = 0.8 (dashed curve), with / e ,iSM — 0.6 (dotted curve), 
with / e ,lSM = 0-4 (dot-dashed curve), and with / e ,lSM =0.2 
(dot-dot-dashed curve) are plotted. The open square is the upper 
limit estimated from the constraints from Ly a damping wing of 
GRB at z = 6.3 (Totani et al. 2006). The filled circles are the up- 
per limits estimated from observations of LAEs at z = 6.6 (Ouchi 
et al. 2010). The filled squares are observational values evaluated 
with the spectra of QSOs at z ~ 6 (Fan et al. 2006). 



cape fractions of HI, Hel, and Hell ionising photons from 
the massive halo shown in section [3.2. H The computed es- 
cape fractions are 0.31, 0.35, and 0.13 for Hel, Hel, and Hell 
ionising photons, respectively. However we still have many 
unknowns, e.g., the sufficient mass resolution to evaluate 
the escape fraction, a typical value of the escape fraction 
of a halo with a given mass, the dependence of the escape 
fraction on halo mass and/or redshift, and so on. We will 
do a statistical study on the escape fraction to clarify such 
unknowns in the future. 



4.3 Optical Depth for Thomson Scattering 

The measurement of the optical depth for the Thomson 
scattering of CMB photons by electrons released during the 
reionisation is frequently applied for evaluating the reionisa- 
tion epoch. The WMAP 7year have shown that the optical 
depth is r e = 0.087 =b 0.014 , which corresponds a reioni- 
sation at z r — 10.4 ±1.2 under the assumption of instan- 
taneous reionisation. In this section, we present the optical 
depth computed from the simulated reionisation history. It is 
likely invalid to compare the optical depth computed from 
our results with the observational value, since the volume 
size of the simulations is insufficient. However comparing 
the optical depths computed from different runs is useful to 




rN256L5RT 
rN128L5RT 
rN128L5RT_fe02— 
rN128L5RT nfsh 



2 4 6 8 10 12 14 16 18 
redshift z 

Figure 13. Optical depths for the Thomson scattering. The solid, 
dashed, dot-dashed, and dot-dot-dashed, curves respectively show 
the optical depths evaluated from the rN256L5RT, rN128L5RT, 
rN128L5-fe02, and rN128L5RT-nfsh runs. 



understand the simulation results. We evaluated the optical 
depths in our simulations by 



: ca t 



f 



n e (z)^dz, 
az 



(13) 



where c, at, and n e (z) are the speed of light, the corre- 
sponding cross-section, and the electron number density 
at a given redshift, respectively. Since our simulations are 
stopped at z = 6.0, we cannot evaluate electron fractions 
at z < 6.0 with the simulations. Therefore, to compute the 
optical depth, we assumed that the H and He in the IGM 
at z < 6.0 is fully ionised. Note that we probably overes- 
timate the optical depth at z < 6.0 with the assumption, 
since Hell reionization was expected to take place at z ~ 3. 
In Fig. 1131 we show the optical depths computed in sev- 
eral runs , and compare them with the WMAP result. It is 
obvious that the optical depths of all runs never reach the 
observed value. Since the first UV sources would appear at 
higher redshift as the mass resolution increases, we simply 
speculate that the total optical depth is higher in a higher 
resolution run. As we can see in Fig. 1131 the above specula- 
tion is indeed correct, but the increment of the optical depth 
is small. This originates in the fact that the ionisation his- 
tory in rN256L5RT run is delayed by UV feedback during 
the later stage (see Fig. [3}. 

As we can see immediately from Fig. 1131 any the op- 
tical depths estimated from our results are not consis- 
tent with observed value. However, no conclusions can be 
drawn from this, since our simulation volume is not suf- 
ficient to directly compare the optical depth. In addition 
our simulations didn't include Pop III stars which are ex- 
pected to have a large contribution to the early stage 
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of the reionsati on (|Ciardi et al.l 120031 : ISokasian et al.l [2004; 
lAhn et alJl20l3 ). 

4.4 Chemical Feedback 

In this paper, chemical feedback by SN explosions was not 
included. The metal enrichment likely plays an important 
role. The increase of metal directly leads to an enhancement 
of the cooling rate owing to additional cooling processes. 
Moreover, H2 molecules are efficiently produced on the sur- 
face of dust grains. This also leads to an enhancement of 
cooling rate below 10 4 K. Not only these effects on cooling 
rates, but also the absorption of UV photons by grains are 
closely related to the cosmic star formation and reionisation 
history. For instance a direct consequence is that the number 
of ionising photons escaping from halos decreases. The ab- 
sorptions of UV photons by grains would alleviates UV feed- 
back, but push the matter outward by the radiation force. 
Hence the net effect of metal enrichment on the star forma- 
tion and reionisation history is still uncertain. Moreover, it 
is expected that the metal enrichment process regulates the 
transition from Pop III to Pop II star formation mode. We 
will study the complicated role of the metal enrichment in 
a forthcoming paper (Paper II). 

4.5 Resolution Effects 

In this section, we discuss how numerical resolution affects 
our results. As we can see from the lower panel of Fig. O 
our simulation result have not converged; The cosmic star 
formation rates at the later stage in the high resolution 
rN256L5RT run is lower than those in t he low resolution 
rN128L5RT run. ISusa fc Umemural (|2004h have shown that 
a considerable fraction of the gas in a small galaxy is evapo- 
rated by UVB, but the star formation in high density regions 
is still possible owing to the self-shielding effect. Since our 
simulations cannot resolve some small clouds where the self- 
shielding effect originally works, galaxies in the low resolu- 
tion run are thought to be more sensitive to the UV feedback 
than those in the high resolution run. Due to this type of 
resolution effect, we might underestimate the baryon mass 
fractions and the star formation rates especially in low mass 
galaxies. However this resolution effect cannot explain why 
the star formation rates at z ^ 9 in the high resolution run 
is lower than those in the low resolution run. 

There is another type of resolution effects in our simu- 
lations. As we increase the mass resolution, the epoch when 
the first source forms is shifted towards higher redshifts, 
since our numerical resolution is not enough to resolve Pop 
III halos which are the first collapsed objects in the Uni- 
verse. This fact indicates that the feedback processes start 
to work from an earlier epoch in the higher resolution run. 
Since a massive galaxy forms via coalescences of less mas- 
sive galaxies which are already damaged by previous feed- 
back episodes, the star formation in the high resolution run 
is likely delayed compared to that in the low resolution run. 
This trend caused by the latter type of resolution effect is 
basically consistent with the trend shown in Fig(3] There- 
fore it is very important to increase the numerical resolution 
up to resolving the t he scale of Pop III halo M ~ 10 5_6 M© 
([Yoshida et al.l l2003h , not only for understanding the con- 



tribution of Pop III halos to the evolution of the Universe, 
but also for removing this type of resolution effect. 



5 CONCLUSIONS 

We explored the impact of UV feedback on reionisation his- 
tory using RHD simulations in which star formation reg- 
ulated by the internal and external UV feedback in each 
galaxy is appropriately considered. As a result we found 
that the cosmic star formation rate is considerably decreased 
due to the UV feedback. In halos with M < 1O 9 M , the 
star formation is mainly suppressed by photo-evaporation 
and photo-dissociation. Interestingly, more massive halos 
also suffers from a strong suppression of star formation, al- 
though they hardly lose their baryon components by photo- 
evaporation. For massive halos, the suppression by photo- 
dissociation is not very effective, but the high density gas 
clumps are significantly damaged by the internal UV and 
SN feedback. In addition, mild photo-ionisation in the mas- 
sive halo turns out to keep the gas hot, and consequently H2 
molecules are destroyed via collisional dissociation. We also 
found that the degree of the suppression is larger in higher 
resolution run, because the UV feedback works since an ear- 
lier stage, and small galaxies that can be resolved only in the 
high resolution run are more sensitive to the UV feedback. 
Therefore the reionisation process turns out to be delayed 
at z ^ 9 in the high resolution run, compared to that in the 
low resolution run. This trend is indeed a kind of the reso- 
lution effect, but this fact is indicating that resolving Pop 
III halos is essential to compute the correct star formation 
and reionisation histories. 

We also find that the simulated cosmic star formation 
rate densities and HI fractions at z ^ 7 — 8 are well consis- 
tent with those inferred from observations, even though our 
simulation volume is too small to be a representative patch 
of the Universe. 
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